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PROJECT  OVERVIEW  -  HEDM 
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Incident  Fluence  (J/cm2) 
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Five-level  model  for  nonlinear  absorption 
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THEORETICAL  METHODS 
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a)  “Self-consistent  field”  (SCF):  reasonably  good  geometries 

b)  “Electron  correlation”:  post-SCF  correction,  required  for  reliable 

energetics  (e.g.,  barriers). 


THEORETICAL  METHODS 


1p3\'x 


mm 


LL 

oo  a  a  a. 

m  -U  T3  13 

^  o  o  o 


r  Q,  CL  Q.  SI 

>  T3  T3  TJ  T3  . 

CD  O  O  O  O 


i  a.  a  a.  Ql 

I  -D  -O  13  .  T3  T3 

D  o  o  o  o  o 


Q.  Q. 
T3  T3 
O  O 


— -  5E  a  a  a  a  Ql 

C)  T3  "D  T3  "O  T3 

3-  Q£  O  O  O  O  O 


2. 

-o  ‘ 

O 


O.Q.Q.Q.Q.QJQ.  Q- 
(/)0£000000000 


c 

a  a 

£5 

^  3 

o 


c  c  c 

o  .5  .5 

H5  w  </)  *-» 

rc  w  </)  c 

o  a>  >j  0 

(5  I  I  P? 


O  W 
>  ’4=  E 
O)  >  0 


^  i*. 

0  0^Q)  +J  ^ 

T"  -r  O)  ■—  .  C  D) 
■*-  -i-  i_  X3  >(i)  s- 

O  o  ©  0  Q>  0 

r-  ■_=  C  &-  b:  T3  C 


^5^0O)|00 


c  Q.  d  I 


LU  <  Z  <  S  IE 


_  _  LL 

O  O  Q 


DFT  gradient  cdp  cdp 

c  =  conventional  disk  storage  of  AO  integrals 
d  =  direct  evaluation  of  AO  integrals 
p  =  runs  in  parallel 
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Speedup  relative  to  32  CPUs 


THEORETICAL  METHODS 


classical  dynamics  done  on  collection  of  quasiparticles  >  natural, 
efficient  parallelism. 
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Number  of  Nodes 


Input>rocket ,  chamber,  20.410000,  exhaust,  0.010000 

Rocket  specific  impulse  calculation: 

The  chamber  pressure  =  20.41  atm 
The  exhaust  pressure  =  0.01  atm 

The  initial  equation  error  was  huge:  24728.149173 
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The  Chamber  State: 


Reference  state  =  reactants 
H (R)  =  H-1773 . 07 ,  E(R)  =  E-1773.05,  S (R)  =  S-  0.00 


P 

(ATM) 

20.4 


V 

(CC/GM) 
909 . 6888 


T 

(K) 

6436.6 


H  (R) 

(CAL/GM) 

0.00 


E  (R)  S  (R)  VGS 

(CAL/GM)  (CAL/K/GM)  (CC/GM) 
-449.64  2.389  909.6888 


Product  concentrations 

Name  (mol/kg)  (mol  gas/mol  explosive) 
n2  Gas  2.112e+001  5.410e+000 

CO  Gas  7 . 734e+000  1.981e+000 

no  Gas  4.609e+000  1.180e+000 

o2  Gas  1 . 561e+00  0  3.998e-001 

co2  Gas  7 . 513e-002  1.924e-002 

no2  Gas  1.418e-003  3.632e-004 

*c  solid  0 . 000e+000  0.000e+000 

Total  Gas  3.510e+001  8.990e+000 

Total  Cond.  0.000e+000  0.000e+000 


The  Exhaust  State: 


Reference  state  =  reactants 
=  H-1773. 07,  E (R)  =  E-1773.05,  S (R) 


=  S-  0.00 


p  V  T 

(ATM)  (CC/GM)  (K) 

0.0  565350.1813  2148.9 


H (R)  E (R)  S(R)  VGS 

(CAL/GM)  (CAL/GM)  (CAL/K/GM)  (CC/GM) 
-1863.53  -2000.43  2.389565350.1813 


Product  concentrations 

Name  (mol/kg)  (mol  gas/mol  explosive) 
n2  Gas  2.337e+001  5.984e+000 

co2  Gas  6 . 157e+000  1.577e+000 

co  Gas  1 . 652e+000  4.232e-001 

o2  Gas  7 . 645e-001  1.958e-001 

no  Gas  1.234e-001  3.160e-002 

no2  Gas  5.283e-006  1.353e-006 

*c  solid  0 . 000e+000  0.000e+000 

Total  Gas  3.206e+001  8.212e+000 

Total  Cond.  0.000e+000  0.000e+000 


The  specific  impulse  = 


402.54  seconds 


RESULTS  -  HEDM 


Computational  requirements:  4500  MW  memory,  -350,000  node-hrs,  ERDC  T3E 


Identifying  a  Completely  New  Molecule: 
Comparison  of  Calculated  and  Measured  Spectra 


tCCSD(T)/6-311+G(2c/)  results 


Identifying  Precursors  for  New  Polynitrogens 
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Computational  requirements:  —10,000  CPU-hours,  1200  MW  on  IBM  SP/P3  at  ASC 


Other  Potential  Polynitrogen  Precursors 

Being  Investigated 


Atom-Doped  Solid  Hydrogen 
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Computational  requirements:  -50,000  CPU-hrs,  MHPCC  IBM  SP 


Atom-Doped  Solid  Hydrogen 


Atom-Doped  Solid  Hydrogen 
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Atom-Doped  Solid  Hydrogen 


Kriimrin.c,  J.R.,  Jang,  S.,  Alexander,  M. PI.,  and  Voth,  Mirijanian,  D.  1 Alexander, 

G.A.:  .1.  Cliein.  Phys.  113  (2000)  9079  published  ;  :  f  Vy  : 


Ring  Expansion:  RSi(OH)3  +  D3  -»  D4  +  H20 

.  Condensation  of  D3,  D4  to  POSS  (in  progress) 

2Da  ->  To  +  4I-LO 
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RESULTS  -  POSS 


Reactive  Complex _ j _ Reactive  Complex 


RESULTS  -  POSS 
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